Spin transistor action from Onsager reciprocity and SU(2) gauge theory 
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We construct a local gauge transformation to show how, in confined systems, a generic, weak non- 
homogeneous 5(7(2) spin-orbit Hamiltonian reduces to two U(l) Hamiltonians for spinless fermions 
at opposite magnetic fields, to leading order in the spin-orbit strength. Using an Onsager relation, 
we further show how the resulting spin conductance vanishes in a two-terminal setup, and how 
it is turned on by either weakly breaking time-reversal symmetry or opening additional transport 
terminals. We numerically check our theory for mesoscopic cavities as well as Aharonov-Bohm rings. 

PACS numbers: 72.25.Dc, 73.23.-b, 85.75.-d 



Transistor action is often based on symmetries. To 
switch on and off a field effect transistor, an exter- 
nal gate turns a three-dimensional insulator into a two- 
dimensional metal and back. Compared to the off-state, 
the on-state has thus reduced dimensionality and sym- 
metry. The relevance of symmetries in transistor ac- 
tion is even more pronounced in some recently proposed 
spin-based transistors, whose action follows directly from 
the breaking of spin rotational symmetry, by tuning 
spin-orbit interaction (SOI) around a special symmetry 
point pi, where the SOI field reduces to two identical 
C/(l) fields with opposite coupling constants 0]. 

In this manuscript, we propose a new class of spin tran- 
sistors whose action is based on an Onsager reciprocity 
relation. We show that in confined quantum coher- 
ent systems with spatially inhomogencous SOI (Rashba, 
Dresselhaus or impurity SOI, or a combination of the 
three), an appropriate SU(2) gauge transformation al- 
lows to express the spin conductance G[- — between 
two terminals labeled i and j through the charge magne- 
toconductancc G VJ (3) as G\- - G^ = G tJ (3) - G ij (- r B). 
This holds to leading order in the ratio L/£ so <C 1 of 
the system size L and the spin-orbit (precession) length 
4o- The pseudo magnetic field 3 arises from the gauge- 
transformed SOI. Current conservation together with the 



Onsager relation G l] {'S>) = G ji (-'B) 0,1] then forces 
G?t- — G\- — to leading order for a two-terminal setup. 
This is the off state of our transistor. The on state is ob- 
tained by cither opening additional terminals or breaking 
time- reversal symmetry with a true magnetic field Bq, in 
which case G\- - G\- = G„-(S + 3) - G y (B - 3) ¥= 0, 
even in a two-terminal setup. Our Onsager spin transis- 
tor can thus be controlled either electrically or magneti- 
cally. In both instances, this turns on a spin conductance 



-4 



with an on/off ratio cx (£ BO /L) 1. 
The mechanism works in diffusive as well as ballistic sys- 
tems, and is more pronounced in regular systems with 
few channels. 

Aleiner and Falko constructed a gauge transformation 



to show that, in confined systems with L/£ so <C 1, a 
homogeneous fc-linear SOI has a much weaker effect on 
charge transport than the naive expectation cx £~^ [H. 
Brouwer and collaborators later argued that terms cx £~J 
in the charge conductance survive the gauge transfor- 
mation for SOI with spatially varying strength The 
relevance of the pseudo magnetic field for a specific meso- 
scopic system with inhomogenous SOI was noticed in 
Ref. 0]. It is however not clear how much of the gauge 
arguments of Refs. JB|, @, Q carry over to spin transport 
in generic systems [9( . Below we show that gauge trans- 
formations result in different symmetries for charge and 
for spin transport [Icj ]. 

Our starting point is a two-dimensional Hamiltonian 
for electrons with SOI, which we write as (H = 1) 



H 



--D^ + V( X ) 



(1) 



Here, V^x) is a spin-diagonal potential and the covariant 
derivative = — (ik so /2)a a A 1 ^ contains the SOI via 
the SU(2) gauge field a a A^, with the Pauli matrix a a . 
Prom here on, Latin indices are spin indices, while Greek 
letters denote spatial indices. The SOI constant fc so de- 
termines the spin-orbit length as £ so = ttI^soI^ 1 - We 
consider a gauge transformation — > 0' = UOU^ 1 with 
U = cxp(icr a A Q /2) ~ 1 + icr a K a /2, and search for a A a 
that reduces the leading order, fc so -linear part of the SOI 
to a spin-diagonal U(l) x <j z structure. We use the well- 
known decomposition (e^ v is the totally antisymmetric 
tensor of order two) for each spin component 



4" 



-{d^ + ^udAl , (2) 

with (j) a given by V 2 a = e^d^A". In particular, 4> a 
is necessarily nonzero for spatially varying SOI. It is 
straightforward to see that the choice A a = fc so x a gauges 
away the gradient part of the SU{2) vector potential to 
linear order in fc so , 
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(3) 
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Note that O(fc so ) corrections in A a u lead to 0(k^ o ) correc- After the global spin rotation er-f -^a z , Eq. JTJ) becomes 



tions in the Hamiltonian. If the SOI strength is spatiall 
constant, <fi a = and one recovers the result of Ref. 
that all O(fc so )-terms are gauged away. 

We next want to extract the leading order, linear in 
fc so X a <C 1 dependence of transport properties such as 
conductances, and thus use 

0' = UOU- 1 = + ik so [a a X a , 0]/2 . (4) 
In particular we have — >• D'^ with 

a -> a =a + k so e x Wcr . (b) 

To calculate spin conductances we need to gauge trans- 
form the operator for spin current through a cross-section 
Cj in terminal j, I" — J c da {n Q ■ j(a),a a }. We obtain 



/ da [{n a ■ j'(a), a a + k so e abc X b (a)a c }} 



(6) 



where j" is the "naive" spin current of the transformed 
Hamiltonian, not accounting for the rotation ([5]) of the 
spin axes. We further need the Heisenberg picture oper- 
ators If (t) = Q lHt l^ e - %m which transform as 



(I')°(t) = ,J*(t) + k so 5J« (t) + 0(k 2 so ) 



(7) 



Here the subscript means that the time-evolution is 
through the fc so = Hamiltonian. 

Linear response relates chemical potentials in exter- 
nal reservoirs and currents in the leads via the spin- 
conductance matrix as If = J2j GfjUj/e. It is some- 
how tedious, though straightforward to show that, to lin- 
ear order in fc so , the gauge transformation gives Gfj — > 
(G^) a , with the conductance matrix (Gy)° evaluated 
in the same way as Gfj but with the spin current op- 
erators Jj of the transformed Hamiltonian in Eq. (JTJ). 
Thus, to leading order in k so , infinitesimal nonabelian 
gauge transformations preserve the form of the spin con- 
ductance. Note that global gauge transformations (i.e. 
global spin rotations), whether infinitesimal or finite, are 
easy to introduce via the corresponding rotation matrix 
ft as Gfj = ft Qb (G^) b . All global or local spin gauge 
transformations leave the potential V(x) invariant. 

We are now equipped to use the gauge transformation 
to explore the spin conductance. We first focus on the 
exactly solvable case of a Rashba SOI [ll| with a spatially 
varying strength a(x) = fc so a (x-f), with a dimensionlcss 
function a, whose gradient always points in the direction 
of the unit vector f . One has — — 2a(x-f )e aAl , </> a (x) = 
ip(x)f a , and Eq. (j3|) gives 



H 



h{a) 
h(-a) 



(9a) 



D' = <9 M + -k &0 e^ v d v ip(-x.)cr-i. 



(8b) 



h(a) = --^[V + jfc s „a(x)f+K(x) . (9b) 

Thus to linear order in fc so , the Hamiltonian is mapped 
onto a block spin Hamiltonian where the opposite spins 
feel opposite, purely orbital pseudo magnetic fields 23 = 
(V x a) z generated by the U(l) vector potential a M = 
^e/i!/<9„</?(x). We obtain S(x) = fe so f • Va. Transforming 
back to the original gauge, the spin conductance is ob- 
tained as G% = [Gij(S) - Gy (-S)]/ a + 0(fc s 2 o ). In this 
simple example, the spin conductance is thus the differ- 
ence of two charge conductances Gij at opposite pseudo 
magnetic fields. For generally varying SOI, one cannot 
choose a spin quantization axis as before. Thus we need 
to define one pseudo-magnetic field per spin polarization, 
i.e. we define S a = d x Ay — dyA® as the magnitude of 
a pseudo magnetic field (pointing always in z-direction) 
that arises solely from the a component of 4> a . To linear 
order in k so L, the superposition principle gives the spin 
conductance along axis a as solely due to the component 
of a , Gfj = G, 3 (S„) - Gy(-® a ). The same argument 
gives the leading-order spin conductance in the presence 
of an externally applied (i.e. true) magnetic field Bo as 

G%{B ) = G l3 (B + S„) - G^ (Bo - S„) + 0(fc s 2 o ). (10) 

This is our main result. It expresses the spin conductance 
of the original dot with SOI in terms of charge conduc- 
tances of the dot without SOI, but with effective mag- 
netic fields Bo ± 25 a arising from the true applied field, 
Bo, and the pseudo field, 23 a , generated by the gauge 
transformation and the SOI. 

The key observation is then that the reciprocity rela- 
tion Gij(B) = Gji(-B) [H, together with gauge invari- 
ance, z^2jGij(B) = 0, imply that the spin conductance 
(fTU|) vanishes to order O(fc so ) in two-terminal geometries 
in the absence of external magnetic field, since only then 
G? ii (-S ) = G ii (S a ) = G ij -(S a ). On the contrary, G% is 
linear in fc so , i.e. much larger, when an external magnetic 
field is applied or when one (or more) additional termi- 
nals are open. Thus, multi-terminal spin conductances 
linearly depend on k ao L, whereas two-terminal local con- 
ductances are quadratic or higher order in k so L. These 
restrictions imply that any coherent conductor with spa- 
tially varying SOI can be operated as a spin transistor, 
whose action is controlled by either opening an extra ter- 
minal or applying an external magnetic field. This is the 
fundamental mechanism on which the Onsager spin tran- 
sistor we propose is based. 

We numerically confirm these results by computing 12 1 
the charge and spin conductances for two- and three- 
terminal mesoscopic cavities and rings (sketched in the 
inset of Figs. Q}[3|). We first assume a Rashba SOI 
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Figure 1: Comparison of the spin conductance G V RL = 
G RL — G RL with the difference in the magnetoconductance, 
Eq. (|f 01) . for transport (from left to right lead) through the 
three-terminal ballistic quantum dot (see inset) with linear 
size L, leads of width W and spatially varying SOI a(x) = 
fcoQ(x) = k so (y/L) (i.e. "B = d y a = k so /L) with k BO L « 0.3. 



with constant gradient over the whole conductor, a(x) = 
(y/L)k so , and check the prediction (|10[) that the spin 
conductance can be expressed in terms of the charge 
conductance of the transformed system without SOI but 
with a magnetic field 23. In Fig. [TJ the spin conduc- 
tance G V RL = G RL — G RL (from now on the y-axis is the 
spin quantization axis) in the absence of magnetic field 
is compared to the difference of the charge conductance, 
Grl{'S')— Grl(— 23) in the absence of SOI, but with mag- 
netic field 23 = d y a. Both quantities exhibit precisely the 
same mesoscopic conductance fluctuations as a function 
of Fermi momentum, as predicted by Eq. (|10p . We found 
that this level of agreement holds up to k so L sa 1, be- 
yond which terms quadratic and higher order in k so are 
no longer subdominant. 

For weak magnetic fields (with an associated cyclotron 
radius larger than L), Gij(B) — Gij(-B) is predomi- 
nantly given by quantum coherent contributions only. 
They give rise, on top of the mesoscopic fluctuations dis- 
played in Fig. [TJ to a shift 8G in the (energy) averaged 
conductance, known as weak localization correction. In 
the presence of a magnetic field, SG exhibits a damp- 
ing that is Lorentzian-like, 5G(B) = <5G(0)/(1 + £B 2 ), 
for chaotic ballistic cavities [HI with (5G(0) ~ (l/4)e 2 //i 
and £ proportional to the dwell time in the cavity. Ac- 
cording to the prediction (|10[) for the two-terminal case, 
the presence of an external magnetic field Bq leads to a 
finite spin conductance G y {B ) = G(B Q +'B)-G(B -'B), 
with 23 = d y a. Then its energy average is 



(G V (B )) 



SG(0) 



SG(0) 



l + £(B + 23)2 1 + £(B -S)2 



(11) 




This line of reasoning is confirmed in Fig. EJa) where 
numerically calculated spin conductances (symbols) for 
the chaotic cavity with linearly varying SOI are compared 
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Figure 2: Spin conductances for two-terminal geometries 
as a function of an applied magnetic field Bq. (a) Average 
spin conductance of a chaotic cavity (inset) for four different 
strengths of a linearly varying SOI (same as inset Fig. [TJ) with 
k so L ~ 0.16,0.33,0.67 and 1.0 from bottom to top curve. 
Symbols with statistical errorbars mark numerical results for 
the average spin conductance, full lines depict the theoretical 
prediction {TTJ). The grey dashed line shows predicted spin 
conductance maxima (from Eq. (JTTJl) for varying gauge field, 
(b) corresponding on-to-off ratios (G y (Bo)) / (G y (0)). (c) spin 
resolved conductances G^'($) for an AB ring (inset panel 
(d)) as a function of flux "I> = 7r_R 2 _Bo, showing a shift ±2? due 
to the gauge field 2 = Va arising from SOI a — (y/L)k so 
with k so L = 1. (d) resulting spin conductance G y (&) of a 
single AB ring. Inset panel (a): Sinai-type billard: linear 
size L, stopper disk with radius Ri — L/10, leads of width 
W = L/15 hosting 4 transverse channels. Inset panel (d): AB 
ring: radius R = L/2, width W = L/15 with 4 open channels. 



to the prediction (fTTJ) (full lines). Figure EJb) shows the 
corresponding on-off ratios (G y (Bo)) / (G y (0)) . 

Alternatively, we consider few-channel regular 
Aharonov-Bohm (AB) rings where fc so -hnear spin cur- 
rents can be turned on by a magnetic flux Q. These 
systems exhibit large almost periodic AB conductance 
oscillations instead of the weaker, randomly-looking 
conductance fluctuations. In Fig. [Uc) we present nu- 
merically computed spin resolved conductances G^* 1 (<&) 
as a function of flux $ = irR 2 B n (in units of the flux 
quantum $ = h/e) for an AB ring (inset panel (d)) in 
presence of the same linearly varying SOI as for the cav- 
ity. As expected, the conductance traces for the spin-up 
and -down channels arc shifted against each other by 
±23 = ±d y a. This shift gives rise to a finite i?o-pcriodic 
spin conductance G y = G^ — G^ as displayed in Fig.[2jd). 
At i?o = 0, first order spin conductance is forbidden 
by the Onsager relation. G y vanishes further for fields 
corresponding to <j> , $o/2 and <&o/4, where maxima and 
minima of the usual charge magnetoconductance occur. 
Maxima of G y appear at points where the shifted spin 
resolved G^^ have their minima. This holds for regular, 
or quasi-regular electronic dynamics which requires clean 
AB rings with few-channels. Of particular interest in the 
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Figure 3: Onsager symmetry-based transistor action re- 
sulting from the difference in spin conductance of two- and 
three-terminal mesoscopic rings (insets panel (b)). (a) On- 
off ratio and (b) separate spin conductances G| T and G\ T 
for AB ring in three- and two-terminal mode as a func- 
tion of a spatially nonuniform SOI, a = (y/L)k so . (c) 
double-log plot of same data as in (b) (top and third sym- 
bol sequence) and of corresponding Gg T and Gj T (second 
and fourth sequence) for a more generic nonuniform SOI 
a — k BO cos 2 (2itx/ Li) cos 2 (2ny / L2) ■ 



AB case arc: (i) the magnitude of the spin conductance, 
which exceeds its value in chaotic systems by one to two 
orders of magnitude (compare the vertical axes scales in 
Fig. [Handed)), and (ii) the control one has over the 
spin conductance: Applying an integer or half-integer 
flux quantum gives the off state of our transistor, while 
the on state is recovered at B Q = ±[(<& /4)/(7ri? 2 ) -23]. 
The on/off spin current ratio can be made arbitrarily 
large, as it exactly vanishes in the off state. 

As said above, fc so -linear spin conductances can also be 
turned on by adding an additional terminal. As shown 
in Fig. Elja,b) we find a difference of at least three or- 
ders of magnitude in spin conductance, G\ T vs. Gg T , for 
two- and three-terminal rings. In panel (c) a double log 
representation of the data from (b) reveal the cubic vs. 
linear k so L dependence of G^p (top symbol sequence in 
(c)) and Gg T (third sequence from top) in line with our 
predictions. 

So far we have considered linearly varying SOI. How- 
ever, our theory holds generally and works well also for 
more generic spatial dependence of the SOI. We confirm 
this by calculating G y (k so L) for a ring with SOI a(r) = 
k so cos 2 (2-kx / Li) cos 2 (2iry / L 2 ) with LjL\ = 15, L/L 2 = 6 
giving rise to SOI bumps on scales of the ring width. As 
demonstrated in Fig. [3]Jc) we recover again the linear vs. 
cubic scaling with k so L for the two- and three-terminal 
setting (second and fourth symbol sequence from top), in 
full accordance with our theory. 

We conclude with a few remarks: 

(i) Mesoscopic rings based on InAs [3| or p-doped 
GaAs samples which are known to exhibit large and tun- 
able SOI [l5| are excellent candidates to experimentally 



probe our theory. Inhomogcneous SOI could, e.g., be 
realized through a top gate covering only part of the sys- 
tem. Additionally, a measurement protocol for spin cur- 
rents based on symmetries of charge transport through 
quantum point contacts [l6| could be implemented. 

(ii) Inhomogeneous SOI is also a prerequisite for vari- 
ous specific proposals for spin splitting EH and ana- 
logues of the Stern-Gerlach effect Our theory pro- 
vides a rather general, common footing to interpret them. 
For instance, the Stern-Gerlach based spin separation, 
usually explained in terms of a Zeeman coupling in a 
non-uniform (in-plane) magnetic field (associated with 
Rashba SOI), finds its explanation in the opposite bend- 
ing of electron paths owing to the Lorentz force associ- 
ated with our gauge field ±23. 

(iii) Another gauge transformation, dual to ours, al- 
lows to transform a nonuniform Zeeman term into two de- 
coupled components with an additional gauge field [2oT ]. 

(iv) While the spin conductance fluctuations are sim- 
ilar in a (phase coherent) diffusive system, its classical 
magnetoconductance has a linear in magnetic field con- 
tribution originating from the classical Hall effect. Thus 
in a diffusive system with inhomogeneous SOI, we expect 
a spin conductance with a nonzero average value propor- 
tional to the classical Hall conductance. This spin con- 



ductance can be estimated [21( as (G a ) 



(e 2 /h)(k so £) 

where I is the mean free path. We stress that (G a ) is 
based on a classical effect in that it is robust against ef- 
fects such as dephasing and temperature broadening. 
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